Dirac lattice 



Corneliu Sochichiu* 
University College & Physics Dept., 
Sungkyunkwan University, Suwon 440-746, KOREA 

November 13, 2012 



Abstract 

We study the emergence of Dirac fermionic field in the low energy 
description of non-relativistic dynamical models on graphs admitting 
continuum limit. The Dirac fermionic field appears as the effective 
field describing the excitations above point-like Fermi surface. To- 
gether with the Dirac fermionic field an effective space-time metric is 
also emerging. We analyze the conditions for such Fermi points to 
appear in general, paying special attention to the cases of two and 
three spacial dimensions. 



1 Introduction 

It is firmly established, that we leave in a relativistic world. There was no 
slightest deviation from the Lorentz symmetry observed so far. The events 
around OPERA experiment recently confirmed the solid status of this fact [T]. 
This means, that any reasonable physical theory, at high enough speeds, 
should be relativistic, must respect Lorentz symmetry, or more generally, 
the Poincare symmetry. But Poincare or Lorentz symmetry is not a tame 
subject, either group is non-compact, which is a source of multiple problems 
in quantum theory. 

These problems lead some authors to question the fundamental role of the 
relativistic field models, as well as, the exact nature of the Lorentz symmetry 
|2). There was a number of proposals for models with no exact Lorentz 
symmetry, but which are supposed to flow to a relativistic system in the low 
energy limit (3||4]. The general idea is that the observable relativistic high 
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energy behavior is only an approximation to a more complicated 'super high 
energy' system with no Lorentz symmetry, with remarkably small corrections 
at the present day available energies. 

The Lorentz symmetry is known to naturally appear in the low energy 
description of some Fermi systems. Thus, Tomonaga-Luttinger hquid in 
one dimension (5|[6], has an effective description in terms of an interacting 
1 + 1 dimensional fermion. In fact, Dirac fermion generically emerges in the 
case one-dimensional lattices. In two spacial dimensions, for long time it has 
been known, that the low energy description of the electronic wave function of 
graphene, one atom thick layer of graphite, is given by a Dirac fermion in 2 + 1 
dimensions (?] (see [sj for a review of properties of graphene). The general 
conditions of emergence of Dirac fermions from two dimensional lattices were 
extensively studied in |9]. The low energy theory of superfiuid ^He is also 



known to be related to Dirac fermion 10 . Some three dimensional substances 



exhibit (pseudo)relativistic Fermion behavior as well 11,12 . Thus, in 13 



the emergence of Weyl particle dispersion relations was studied in connection 
with three-dimensional crystallographic groups. 

In this work we address the problem of finding such microscopic systems, 
which in the low energy limit flow effectively to the Dirac fermion system. 
The starting point is a non-relativistic system of particles hopping on a graph. 
These particles are obeying a Pauli exclusion principle. The exclusion prin- 
ciple is essential here, since it is the dynamics near the Fermi point, which 
brings us to the Dirac particle in the low energy limit. We present here a 
simple argument, showing that the spectrum of fiuctuations around a non- 
degenerate linear Fermi point is necessarily described by Dirac particles. 
The only assumptions we make, are the analyticity of the energy spectrum 
near the Fermi point and irreducibility in the continuum limit. The last as- 
sumption is equivalent to emergence of continuous rotational symmetry. The 
fact, that the spectrum of fiuctuations near a Fermi point is described by a 
(pseudo)relativistic fermion, was claimed earlier [l^, based on the complete- 



ness of the ABS construction 14 , but a direct proof was not known to us. 
On the other hand our proof presented here, is direct and rather simple. As 
a consequence, our task becomes to look for a system with a non-degenerate 
spectrum around Fermi points. We can express this as a condition imposed 
on the adjacency matrix describing the graph, in addition to the condition 
that the graph admits a continuum description in the low energy limit. As 
a result, we present a general solution in two and three spacial dimensions, 
as well as a solution in arbitrary dimension, which we believe is general, 
although we could not find the proof of that fact in this work. 

The plan of the remainder of this paper is as foUows. In the next section, 
we set up the problem and introduce the notations. Then, in section |3| we 
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show that non-degenerate fluctuations around a Fermi point in the linear 
order in momentum are described by a Dirac fermion or a collection of Dirac 
fermions. In the last case, the energy spectrum of fluctuations belongs to 
a product of inequivalent Clifford algebra representations. In this case, the 
rotational symmetry in the continuum limit is generally broken. Restricting 
to a single copy of such representations, we restore the rotational symmetry. 
There we also show, that an effective (Minkowski) geometry emerges in this 
limit. The conditions, that we stay on a Fermi point, are formulated in the 
form of algebraic equations (3.16), whose solution is considered in section|4| 
There we give a general solution for the expansion coefficients in the cases 
of graphs fiowing to two and three dimensional spaces. We present also a 
solution for the general case of arbitrary dimension. Finally, in section |5| 
we discuss the obtained results, point to perspectives and generalizations as 
well as further problems. In the appendix, we construct the the holomorphic 
Clifford algebra basis, which is used to expand the adjacency matrix from 
the main part. 



2 The model 

Let us consider a the model of a fermionic particle hopping on a graph. The 
energy of hopping between two sites, say x and y is given by the element 
t<xy> of the adjacency matrix T. Further properties of the graph in the form 
of restrictions imposed on the adjacency matrix we will specify later. So far, 
the most general tight-binding Hamiltonian describing this system is given 
by, 

H = - T ■ SL= ^ t<xy>4S' (2-1) 

<xy> 

where a J. and are, respectively creation and annihilation operator s of a 
fermionic particle at the site x. In what follows we assume that the Fermi 
sea is filling half of the available energy states of this model. 

The hermiticity of the Hamiltonian requires that the the amplitudes sat- 
isfy t^xy> = t*<yx> where the star denotes the complex conjugation. On the 
other hand, the complex phase in the transition amplitude t^xy> can be re- 
lated through Peierls substitution to a background (electro)magnetic field, 
therefore in the absence of such a field it is natural to restrict the model to 
real valued transition amphtudes, t^xy> = t*<xy> — ^<yx>- 

In order to admit a continuum description in the low energy limit, the 
graph should possess at least an approximate translational symmetry. As 
a first approach, here we consider graphs possessing the periodicity of a 
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D-dimensional lattice, leading to the flat space continuum limit. The peri- 
odicity property of the graph is reflected in the foUowing block structure of 
the adjacency matrix T: There is a primitive cell subgraph which consists of 
p vertices, which are periodically reproduced in the form of a D-dimensional 
lattice with lattice vectors i The primitive cell itself is an abstract graph, 
not necessary a part of a physical lattice. The adjacency matrix block cor- 
responding to primitive cell is a p x p matrix M. In addition there are 
transitions between neighbor primitive cells n and n + ?, for i = 1, . . . ,D. 
These transitions are given by D blocks Fj also of size p x p. 

In this work we mostly keep the primitive cell index as an implicit matrix 
index, while the lattice coordinate n will be shown explicitly. 



Taking all above into account, we can rewrite the Hamiltonian (2.1 ), using 
explicitly the block structure of the adjacency matrix, 



n, i 



(2.2) 



The Hamiltonian in the form of (2.2) contains a discrete translational 



symmetry and, therefore, is suitable for Fourier transform, 

""-y, (2^ 



a(K) 



E 



JKn 



:a(K)e 



iK n 



(2.3) 



where the momentum k is integrated over the first Brillouin zone B, defined 
by 

K = Ki^, -Ti <Ki< TT, (2.4) 
where the tilded vectors iA = l D form the dual lattice basis, 



In terms of Fourier transform the Hamiltonian takes the form, 
d^i^ 



H 



(F.e-- +Fte--^) +M a(K) 



(2.5) 



(2.6) 



where the Fourier modes a (K) and a^(K) are a p-dimensional vector and a 
co-vector, respectively. The quantity in the parentheses. 



fe(K) 



E (r 



+ r e 



+ M 



(2.7) 



is called energy matrix. The eigenvalues of this matrix give the energy modes 
of the Hamiltonian. The properties of /i(K) is the main subject of the next 
section. 
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3 Low Energy Limit: Fermi points 



The low energy behavior of the model is determined by the regions of the 
momentum space near the Fermi surface. Therefore, the geometry of the 
Fermi surface is essential for the properties of the model this limit. In a 
generic situation, the Fermi surface is a surface of co-dimension one, i. e. 
a {D — l)-dimensional surface in the D-dimensional space. However, the 
Hamiltonian may contain terms, leading to degeneracy of the Fermi surface 
down to a lower dimensional surface or a gapped state. 



The Atyah-Bott-Sapiro construction 14 , or just symmetry arguments 



imply that the emergence of Dirac fermion in the low energy limit is condi- 



tioned by the presence of a Fermi point 10 . Therefore, let us concentrate 
on finding the conditions of appearance of Fermi points and Dirac fermions. 
Most generally, a Fermi point K* is defined by two conditions, 

det h{K*) = O, det[/i(K* + k)] ^ O, for k ^ 0. (3.1) 

The first equation means that there are zero modes at K*, while the second 
one implies that this is an isolated point. More specifically, we consider 
linearly non-degenerate Fermi points, which imply that for given choice of 
coordinates in the momentum space, the energy matrix is non-degenerate in 
the linear approximation, i. e. 

det[h{K*) + ai{K*)ki] ^ O, for k O, (3.2) 

where «i (K*) = dih(K*). 

In principle, the energy matrix h(K) can be diagonalized, which gives 
the energy bands e„(K), where n = l,...,p. The Fermi point condition 



(3.1) implies that some of these bands must cross zero at the Fermi point 
K*, while other bands may have a gap. In the low energy limit the states 
belonging to the bands having gap will be frozen. Therefore, only gapless 
bands are relevant for this limit. So, one can locally project from the total 
p-dimensional "internal space" to the subspace of gapless modes only. The 
resulting geometry will be generically a fiber bundle, or a K-theory element. 



as described in 15 . Since the contribution of gapped bands is trivial in the 
low energy limit, in this work we assume a minimal case with no gapped 
bands. 

If we expand the energy matrix in the vicinity of the Fermi point, 

h{K* + k) = ai{K*)ki + 0{k^), (3.3) 
where the matrices «j (K*) are given by, 

a,(K*) = i (r,e^^* - rle-'^*) . (3.4) 
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then for small enough nonzero ki, the condition (3.1) becomes, 



h(K*) = O, det[a\K*)ki] ^ 0. (3.5) 

One might suspect, that matrices a^, satisfying the non-degeneracy condition 
(3.5) should be related to Dirac matrices. Inverse is definitely true. In what 
foUows, let us study the properties of these matrices and show that it is 
indeed the case. 



3.1 The properties of matrices ai 

The problem of finding a non-degenerate set of D matrices «j, i = 1, . . . ,D 
can be related to the geometrical problem of finding the maximal number 



of linearly independent vector fields on a sphere 16 . The solution can be 
expressed in terms of Clifford algebras by Atiyah-Bott-Shapiro (ABS) con- 
struction [Til (see also [171 ). 



Still there is a simple way to see, that the second condition (3.5) implies 
that the matrices are generators of Clifford algebra. Consider the eigen- 
value problem for a^ki. It is not difficult to see, that in the vicinity of Fermi 
point the leading term in k should take take the form, 

en{k) = ±\U-k\+0{k'), (3.6) 

where ^„ are some non-degenerate Z)-dimensional square matrices. The ab- 
solute value denotes usual Euclidean norm, 

\y\ = Vy^- 

Restricting to the subspace corresponding to one we have, 

el^{^^-kr + ..., (3.7) 

which means that the leading term of the restriction of the energy matrix to 
the subspace of fixed ^„ must be proportional to the identity matrix (in cell 
index) , 

{a\{K*)hY = {^n-ky = gi^hk,, (3.8) 

where we introduced the n-th subspace "metric tensor" given by the matrix 
9ri = in ' in aud is the restrictiorj^ of matrix a* to the n-subspace. Fur- 



thermore, the equation (3.8) implies that the matrices a*, when restricted to 



n-subspace, should obey Clifford algebra anti-commutation relations, 

a>{ + = 2g^, (3.9) 



^It is a simple exercise to show that the n-subspace is an invariant space for each a% 
therefore such a restriction is well-defined. 
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where gl^ is a positive definite matrix in ij indices. This means that in a 
general basis matrices aj satisfy the anti-commutation relations, 

{a\a^} = 2G'^ (3.10) 

with diagonalizable p x p matrices i, j = 1, . . . ,D with eigenvalues gl^ . 

In the case of distinct gnS, the D-dimensional rotational symmetry re- 
mains broken even in the continuum hmit unless all gl^ are equivalent up to 
a scale factor, gl^ = g^Hf^. The scale factor as well as the metric gij should 
be positive in order to keep the positivity of the energy-square. This resuhs 
in the structure of G*-^, 

where T is a x matrix with the eigenvalues tn in the diagonal basis. 
This implies that in some bases we can split the matrices «j according to, 

«i = «j ® T, (3.11) 

where ai realize an irreducible representation of the Z)-dimensional Clifford 
algebra, 

{a\a^} = 2gij. (3.12) 

Now let us recall, that in the absence of gapped modes, the dimension of 
aj-matrices is given by p, the number of elements in the primitive cell. Since 
the structure of «j-matrices is a product of a Clifford algebra-valued matrix 
to an arbitrary non-degenerate hermitian matrix T, we can evaluate the 
required number of elements in a primitive cell which can yield a relativistic 
fermion in the low energy limit, using the dimensions of respective Clifford 



algebra representations 17 



Thus, for D-dimensional Weyl fermion one has, 

for some positive integer number A^. Dirac representation implies. 

Minimal case would be when = 1. As we will see below, the Fermi 
point in our model is always Dirac point, unless additional constraints are 
imposed. Therefore, the necessary condition for a graph to be Dirac lattice is 
to contain at least 2^^^'^^ elements in the primitive cell. There are not many 
Bravais lattices satisfying this criteria, but detailed analysis of possibilities 
goes beyond the scope of this work. 
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3.2 ClifFord algebra basis and the low energy limit 

As we just established, the Fermi point condition implies that the matrix- 
valued coefficients aj(K*) should be generators of a D- dimensional Chfford 
algebra or a direct sum of D-dimensional Clifford algebras. Furthermore, the 
SO{D) symmetry in the low energy limit requires that it should be copies of 
the same Clifford algebra. 

For the further analysis of the structure of matrices «j (K*) it is convenient 
to expand them in terms of a Standard complex basis of a D'-dimensional 
Clifford algebra, consisting of matrices S/, J = 1, ... D' /2, satisfying the 
foUowing anti-commutation relations, 

{S,, Sj} = {sl, = o, {S,, = 6ui. (s.is) 

In addition, each S/ is a matrix with real entries. It is clear that the even 



dimension D' of the Clifford algebra (3.13) should be equal or greater than 



the lattice dimension D. Therefore, the minimal value of D' is given by, 
-^min — 2[-D/2], where [■] denotes the ceiling function. The complete basis 
of matrices satisfying these conditions is constructed in the Appendix |Aj 

It is not difficult to convince oneself that the matrices T i too, can be 
expanded in terms of basic matrices S/ with real coefficients, i. e., 

T, = Tu^i, rl = Tu^\. (3.14) 

Next we can observe, that the consistency of the Fermi point condition 



(3.5) implies that the intracell matrix M should be expandable in terms of 
hermitian parts of S/, 

M = m/(S/ + Sl), (3.15) 

where mj are real. 

In terms of the Sj-basis the Fermi point equation takes the form, 

hi{K) = 0, I = l,...,D'/2, (3.16) 

where, 

hj{K) = J2^ue''''+mi. (3.17) 



The coefficients Tn should be regarded as parameters of equations (3.16) 



allowing only point-like solutions for the momentum K. It is clear, that if K* 



is a solution to the Fermi point condition (3.16), then, due to the symmetry 



^We use convention, according to which summation is assumed over the repeated index 

I. 
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with respect to complex conjugation, —K* is a solution too. These two 
Solutions carry opposite topological charges |18|, such that the net charge for 
the pair vanishes. Therefore, such a pair represents a minimal set of Fermi 
points. As a degenerate case one can have a single point K* = O with zero 
charge, but such a configuration may suffer from topological instability. 

We postpone the further discussion on the problem of finding appropriate 
Fj/ and m/, till the Section |4| So far let's assume that the equation (3.16) 
admits a minimal solution with a single pair of Fermi points, which we will 
denote ±K*. 

In the vicinity of either Fermi point ±K* the matrices Ui take the form, 

«,(±K*) = iVu (^cosir;(S, - S|) ± sini^;i(S7 + , (3.18) 

where k is the offset momentum in the vicinity of ±K*, i.e. K = ±K* + k. 
The energy matrix in these vicinities takes the form, 

/i±(k) = /i(±K* + k) = a,(±K*)fci + 0{k'^). (3.19) 

Introducing the subspace associated with the sign ± of the Fermi point, 
we can re-write the matrices (3.18) in the foUowing form. 

Oli = Tij (cosK*i{J:j - Sj) ® I - smK*{J:i + Sj) ® ^3) ■ (3.20) 

Thus, matrices «j are linear combinations of Hermitian matrices /3a, a = 
1, . . . , D', defined as, 

/327-1 = -(S, + S|)®cr3, /32/ = i(S/-Sl)®I, I = 1, . . . , D' /2, (3.21) 

&^ = iWa. (3.22) 



I.e. 



where, 

ef-^ = r,,sini^;, ef = r.7cosK;. (3.23) 

The matrices [i a realize a Standard D'-dimensional Clifford algebra basis, 

{/3,,/3fe} = 25„a. (3.24) 

Therefore the coefficients can be regarded as vielbein coefficients for the 
embedding of a D-dimensional plane into MP , where the Clifford algebra 
(3.24) is defined. This embedding produces an induced metric on the D- 

■a 
i 1 



dimensional plane, given in terms vielbeins ^.'^ 



= (3.25) 
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The embedding coefficients and, therefore, the induced metric (3.25) are 



inherited from the adjacency matrix elements t^xy> in the Hamihonian (2.1 ). 



Therefore, even in the case when the graph is a lattice embedded in the Eu- 



chdean space, the 'induced metric' (3.25) does not necessarily coincide with 



the real space metric. They coincide, however, if the transition amphtudes 
on the graph are inversely proportional to the real space distance between 
the nodes, at least in the leading order. This is often the case for the lat- 
tices of physical substances, consisting of single chemical component, like 
e. g. graphene, where this condition is satisfied and the induced geometry 
is indeed equivalent to the real space geometry, up to a scaling factor given 



by the Fermi speed. In more general case, the induced metric (3.25) can be 
interpreted as an anisotropic or tensor valued 'Fermi velocity'. 

To obtain the Standard Dirac fermion action, let us take the foUowing 
steps. First, introduce an 'orthogonal basis' with respect to the induced 



metric (3.25). If we choose the Cartesian basis of the embedding space such 
that basis vectors with numbers a' = 1, . . . ,D span our D-dimensional space, 
while a± = D + 1, . . . , D' belong to the orthogonal completion, then the 
Cartesian momentum in this basis takes the form, 

g'^'=er'fci, a' = l,...,D. (3.26) 



In terms of the Cartesian momentum the low energy Hamiltonian (2.6) 
takes the form, 

h = jI ^^^(g)/3a'g">(g), (3.27) 

where J is the square root of the determinant of the induced metric 



J= y'detll^?,,!!, (3.28) 

and ip{q) and ip'^^g) include only the low energy modes of, respectively, a 
and a^. Matrices /3" belong to the reduction of the D'-dimensional Clifford 



algebra (3.24) to the _D- dimensional subspace. 

Taking into account the commutation relations, we can write down the 
corresponding low energy effective action, 

Sie = J f ^^i^Hgmg) - f dtH{^\i,). (3.29) 



J (27r)^ J 

In the low energy theory one can extend the momentum integration to 
entire momentum space M^. Making the inverse Fourier transform, we arrive 
to real space low energy action, 

5ie = / d^+V (i^^^ - i^^P'^'da'A , (3.30) 
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where da' = d/dx"'' is the partial derivative with the respect to the induced 
'Cartesian' coordinate. Also, in order to absorb the Jacobian factor J, we 
rescaled the field if) — )■ J'^^'^'ip. 

As the next step, let us generate the Dirac matrix 7° from matrices /3" 



given by (3.21), as follows 



/ = i^'/2^i/32---/3^,. (3.31) 
By the construction this matrix anticommutes with all /3", as well as satisfies, 

(7°)2 = -1. (3.32) 
Then, define the Dirac conjugate, 

^ = ij^j^, (3.33) 
as well as the spacial Dirac matrices 7"', 

Once this is done, the low energy effective action takes the 'genuine' 
relativistic Dirac form, 



^ie = -i J d^+ix^7'^a^^. 



(3.35) 



Let us note, that the low energy effective action (3.35) is written in the 
orthonormal frame with respect to the induced metric ( |3.25 ). The Dirac 
matrices in this frame satisfy canonical anti-commutation relations. In the 
case in which the graph is a lattice embedded in real space we can go to the 
real space Cartesian basis using usual basis transformation rules. In general 
case this can lead to an additional Jacobian factor and nontrivial vielbein 
coefficients ^f, resulting in anisotropic Fermi velocity. 



4 Coefficients F^j 

In this section we discuss parameters Fj/, specifying the Fermi point condition 



(3.16). We shall establish the range of appropriate coefficients Fj/, leading to 



point-like solutions for (3.16). The appropriate parameters Fj/ and m/ can 
take values within a range, which can be called moduli space. Each connected 
component of the moduli space corresponds to an equivalence class of physical 
models which can be deformed into each other by continuous modifications 
of the elements of the adjacency matrix. 
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Figure 1: The {D + l)-gon representing the equation ^jTi/e'^' + m/ for 
m/ 7^ 0. The lengths of the sides are given by the angles to the 

horizontal are Ki. In the case of mi = O the {D + l)-gon degenerates to a 
D-gon. (a) A heptagon for the problem in D = 6. (b) A set of two heptagons 
with parallel sites, solving the Fermi point problem can be regarded as finding 
a rigid hinge mechanism. In this case one more independent heptagon is 
needed to make the mechanism rigid. 



The equation (3.16) is homogeneous, therefore the consistent Tn and m/ 
are determined up to a scale factor. Therefore it is convenient to fix the value 
of one non-zero parameter, say m/, for every /. 

For a given /, one can represent graphically the equation (3.16) as a 
{D + l)-gon on the complex plane, consisting of D sides of lengths \Tij\ 
and orientation e'^' as well as one horizontal side of length \mi\ (see Figjl^ 
for an example). So, in the polygon represent at ion, the upper case index 
J = 1, . . . , D' /2 enumerates polygons, while the lower case latin i = 1, . . . , D 
sides inside one polygon. 

The appropriate coefficients Fj/ and mj are those for which all D' /2 poly- 
gons close by unique, up to discrete transformation, combination of angles 
Ki. A continuous freedom in changing the angles K*, would mean that the 
respective Fermi variety is not a point but a higher dimensional structure. 

It is useful to point to a mechanical analogy in terms of arm-and-hinge 
mechanisms, following from the polygon representation of equation (3.16). If 
we normalize the polygons by rescaling their size by mj^, such that they will 
have a common base, we can build a mechanism made of arms of lengths Fj/ 
connected within a polygon by hinges allowing free motion within the two 
dimensional plane. The arms of different polygons having the some number 
i are kept strictly parallel because of common orientation Ki, see Fig. [T]d. 
The "correct" mechanism is one which is rigid, i. e. one in which no angle K^ 
can be modified without breaking at least one arm. 

Let us evaluate the minimal number of independent polygons which are 
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Figure 2: D = 2 situation. There is a unique triangle you can construct 
with given three site lengths. The angles of the triangle are related to the 
momenta in the foUowing ways: a = n — Ki, P = 71 — K2 + K1 and 7 = 71+K2. 

necessary to fix the continuous freedom of Ki. Any polygon, consisting of 
(D+l) edges of fixed lengths can be deformed in the two- dimensional plane by 
{D — 2) independent ways using its vertices as joints. Therefore, one polygon 
condition fixes two out of D angles Ki. To fix all D values uniquely, we need 
D /2 independent polygon constraints if dimension is even, or {D + l)/2 ones 
if the dimension is odd, i. e. the number is D' /2. It is interesting to point 
that this matches perfectly the dimension of the minimal Clifford algebra. 
This match appears both miraculous and beautiful since there is no known 
direct relation between the moduli of complex plane polygons and the Clifford 
algebras. 

In the foUowing subsections we consider first solutions for two simplest, 
but most interesting cases: D = 2 and 3, then we turn to the case of an 
arbitrary dimension. 

4.1 The case: D = 2 

For D = 2, a single equation ( |3.16 ) fixes the Fermi point, 



^ Tiie^^» + mi = 0. (4.1) 

i 

The holomorphic Clifford algebra generator is Si = where 0"+ is the Pauli 
matrix a+ = |(£i + io"2). 

Equation (4.1) implies that two vectors F^ie'''^* form, together with the 



mi-site in the real direction, a triangle on the complex plane, see Fig. |2j The 
lengths of the three sites determine the orientation angles Ki, uniquely, up 
to a reflection. Thus the elementary triangle sine rule reads: 

Fii r^i mi ^^2) 



sin K2 1 I sin Ki \ \ sin(i^'2 — Ki ] 
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Then the choice, 

mi sin K2 mi sin 
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r^i = r '^r ' (4.3) 



sin(i^* -K*)' sm{K; - Kl 



results in a system with two Fermi points located at ±K* = -^{Kl^K2)- 

Let us turn to the low energy fluctuations. The the low energy fluctua- 
tions around the Fermi vacua ±i^* are described by the effective Hamilto- 
nian, 

H,, = j d^k (iji{a+ ■ k)ij+ + ^l(a_ ■ k)ij_^ , (4.4) 

where, 

sin sin / , , ^ 
a±-fc = ±m-^^^-^(-fci + fc>i 

771 

+ .^r.fK'* - K^*^ [-{cosKl smK;)ki + {cosK; smK*)k2] a^. (4.5) 

Introducing the index related to the Fermi point sign ±, the Hamiltonian 
is cast in the form, 

i/ie = j d^k^^a ■ k^, (4.6) 

where 

^ , sin k: sin / , , ^ 

"■^ = "^sin(ir2*-^r)^^ ^ + 

TTl 

+ . _ [-(cos Kl sin K;)ki + (cos sin K^k^] ® I. (4.7) 
sm(^A2 A^j 

The 'Cartesian momenta' and are given by, 

Tfl 

qx = -^J77i 77^ [-(cos Kl sin K*)ki + (cos K^ sin Kl)k2] , (4.8a) 

sin(^A2 A^j 

iiinKl sin /^2 /, , \ 01 n 

■'> = '" Bin(4-A-r) "^'-'^''- 

Introducing Dirac matrices and going back to the real space action as 



discussed in Section |3j we obtain the Dirac fermion action (3.35) in 2 + 1 
dimensions. 

In fact, this general solution f alis in the universality class of graphene. 
Thus, in the case of graphene the Fermi point is located at K* = (27r/3, 47r/3). 
This corresponds to the equilateral triangle: Fn = r2i = mi = t, where t is 



the nearest neighbor transition amplitude (see e. g. [8] and 19 for a review). 
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Figure 3: The two-dimensional Brillouin zone. The gray shaded area is the 
forbidden location for the Fermi point. Since the models with Fermi point 
locations K* and —K* are physically identical, only one triangular zone gives 
distinct graph model. The system becomes degenerate along the boundaries 
between allowed and forbidden zones as well as at the vertices. Due to 
periodicity the left side of the Brillouin zone should be glued to the right 
side and top to the bottom. Because of this the vertices of both triangles 
are, in fact, identified. 



Since in the case of graphene the transition amplitudes are obeying the same 



type of symmetry as the physical lattice, the induced spacial metric (3.25) is 
proportional to the real space metric, 

9ij = vlSij, (4.9) 

where vp = 3ta/2 is the Fermi velocity in graphene, and a is lattice spacing. 

The solution, we constructed in this section, has a set of singular points. 
The location of the Fermi points can not be arbitrary. By inspecting the 
possible orientations of the triangle edges in Figj2j we come to the conclusion, 
that the momenta are constrained to the foUowing regions, 

O < K* < TT, , {-7i<K;<0, , 

^ ' and < ^ ' 4.10 

-71 < k; <-7l + Kl, [7T + Kl <K*<7f. 

These two regions are connected at their corners, taking into account the 
periodicity of the momentum space, see the Fig|3j In fact, since the models 
with Fermi point locations K* and —K* are physically equivalent, indepen- 
dent Solutions are coming from only one of the triangular zones. 
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The solution becomes singular as the Fermi point location approaches an 
edge or a vertex of the allowed zone. Let us analyze these points in more 
details. 

A manifest singularity occurs along the edge for which sm{Ki — K2) = 0. 
If neither of sin K* vanish, the Fermi point is possible only in the case when 
nii = 0. In fact, this solution can be continuously approached from the 



regular one, given by (4.3), by taking the hmit 



mi -^tsm{Kl - k;), sm{Kl - K*) 0. (4.11) 

where t is the scahng parameter. 

Applying the hmit ( |4.11 ) to the expression for the a ( 4.7[ ), we get, 



d ■ k = t{ki — {sin^ K*ai (g) (J3 + sin K* cos K*a2®l}, (4.12) 

where we have chosen the Fermi point to approach the location K* = {K*, K*± 
tt), with K* 7^ O, ±7r. It is clear that the low energy limit of this case describes 
a one-dimensional infinitely degenerate Dirac fermion. The degeneracy of ki- 
netic term is associated with respective degeneracy of induced metric. 

Now, let us switch to the case when the Fermi point approaches one 
of the other two sides of the allowed zone. This corresponds to the limit 
sin K* — )• O for one of i = 1, 2. The parameter m/ remains finite in this case, 
but the triangle in the Fig|2] should flatten to the real axis, which means 
that the other component of the Fermi point momentum should also satisfy 
r j sin K* = 0. This means that unless the Fermi point is touching a vertex, 
the respective parameter Fj should vanish. It is clear that vanishing of any 
parameter Fj means that any dependence on the momentum component fcj 
is automatically excluded in the low energy theory and the model is also 
one- dimensional . 

Now let us look what happens at the corners, i. e. when the sins of both 
Fermi momentum coordinates vanish. In this case one can observe that the 
conjugate Fermi points K* ^ —K*, merge into a single point. Therefore, 
the Fermi point degeneracy, which leads to doubling of components of the 
fermion disappears as well. The solution (4.3) is replaced in this case byj^ 



ri = e, r2 = m-e, (4.13) 
where —00 < ^ < +00. The low energy fluctuations are described by, 

a-k= (^/£1 + (771-0^2)^2. (4.14) 



•^Without loosing generality, we can take K'i = — if| ~ tt, but allow Ti to be both 
positive and negative. 
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This corresponds to one dimensional Dirac fermion model as well. 

Let us note however, that due to fusion of Fermi points in the case of 
the corner location, the protection mechanism due to conservation of the 
topological charge is lost. Therefore, the system is unstable with respect to 
small perturbations, and in the absence of a different mechanism of protection 
will be totally gapped in the case of interactions turned on (see [l8]). 

Let us conclude this subsection by pointing out that once the Fermi point 
location is inside of the allowed zone the model is not much different from 
the tight binding model for graphene (see fs] for a review), and falls into the 
same low energy equivalence class with it. In fact, the only physical Bravais 
lattice, which belongs to this solution is the hexagonal lattice. 



4.2 D = 3 

In three dimensions, the polygon equations (3.16) in general case are given 
by two quadrilaterals, 



3 



^ Fi/e'^» + m/ = O, / = 1,2. (4.15) 



=1 



For non-zero mj, let us normalize the coefficients Tn, by dividing them 
by m/, such that they satisfy the normalized equation, 

3 

J]Fi,e'-^' + l = 0. (4.16) 

i=l 

The original normalization can be restored by redefinition: F j/ — )■ mjTii. 

Observe, that since the edges corresponding to the same direction have the 
same orientations on the complex plane, one can obtain both quadrilaterals 
from a single master triangle, by cutting one angle of this triangle by parallel 
cuts, passing at different distances from the vertex, see the Figjl} This angle 
can be top or bottom, depending on angles and K^. Then the problem 
of finding appropriate Tn, reduces to solving the following equations: 

i. Master triangle equation, 

Fi,e'^^ + F3.e'^^ + 1 = 0. (4.17) 



ii. Cut edge equations, 

Viie'^' + V2ie'^' + Vsie'^' = 0. (4.1^ 
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Figure 4: The quadrilaterals can be obtained from a single master triangle 
with sides 1, T[ and Fg, by cutting the upper angle by the side F27-e'^^. The 
dotted segments correspond to r/i ^. 



We formally wrote two cut edge equations (4.18), but since they are homoge- 



neous, only one of them is independent. As a result, the parameters r]ij can 
differ by a scahng factor independent of space index, r]ij = \ir]i. Therefore, 



we can drop the capital index J in the cut edge equation (4.18). 



Once the master triangle and the cut edge equation are solved, the pa- 
rameters Fj/ can be found from Ti^, and rju = Xi7]i as follows, 



F./ 



1,3, 



21 



(4.19) 



This solution comes in complete analogy with the two-dimensional case. 
Thus, for the master triangle equation we have the solution, 



while for the cut edge equation we can write down the following solution 



(4.20) 



ri2i 



A/, 



XiiAn{Kl-Kl) 
sm{K* - Kl) 



MK* - K*) 
Combining all together, we have, 

m J sin + Xj sm{K2 — 



^21 = -V2I = -A/, 
F37 = miTs^ - r]3i -- 



-m J sin Kl + A/ sin(_ft'^ — K^] 



(4.21) 

(4.22a) 
(4.22b) 
(4.22c) 
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Using the solution (4.22), we can write down the induced metric and 
vielbein coefficients C,i given by equations (3.25) and (3.23), respectively. 
Due to the triangular origin of this solution its properties in the Ki — 
plane mostly mimic that of two- dimensional solution. 

Below we construct the solution based on elementary triangles in an ar- 
bitrary number of dimensions. 



4.3 Holomorphic solutions for arbitrary dimension. 

Inspired by the subtle relations of the polygon equations to Clifford algebras, 
let us find a solution valid for arbitrary dimension. As above, even and odd 
dimensions appear differently: The odd dimensional cases are essentially 
dimensional reductions from higher dimensional even cases. 



Consider the general polygon equation (3.16), 
^r,,e^^'+m, = 0, 1=1,. 



D'/2, 



(4.23) 



and the foUowing Ansatz, 



Ti/ — V'j5i^2i-i + r'/(5j, 



21- 



(4.24) 



As a result of substitution, the polygon equations (4.23) split into D' /2 in- 
dependent triangular equations. 



r,e 



+ r'/e'^^^ + m/ = 0. 



(4.25) 



This Ansatz mimics the canonical form of the rotational matrix, which in an 
appropriate basis is a composition of elementary rotations of two-dimensional 
planes. 

The general solution to equations (4.25) we know already. It is given by 
:ms of nota 

mi sin 



equation (4.3). In terms of notations of this section it spells, 

m/ sin K 21^1 



r'. 



where / 



sin(i^', 
D'/2. 



21-1 



r'; 



21 



^21-1) 



(4.26) 



This solution leads to the following vielbein coefficients (3.23) 
27-1 m/ sin KIj_ ^ sin K 21 



(:2I-1 
?2/ 



sin(i\:* 



2/-1 



K. 



21) 



m J sin K2i_i sin K 21 
■ sin{K*2,_,-K*2,) ■ 



pil 

S2/-1 

p2I _ 
S2/ — 



m/ sin K21 cos 



21-1 



sin{K*2j_,-K*2,) ■' 

mi sin -^'2/-! "^os ^21 
' MK*2j_, - K*2j_,) ■ 



(4.27) 
(4.28) 
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where I = 1, . . . , D' /2. The low energy induced metric has a diagonal 2x2 
block structure, with the block g^j^ given by, 



m, 



9(1) 



X 



sin^ k; 



- sin Kl sin K* cos(ir* - K;] 



sin k; sin K; cos(K* - K;] 
sin^ k; 



. (4.29) 



This solution depends exphcitly on D' /2 + D parameters mj and Ki 



The three-dimensional solution from the Section 4.2, can be obtained from 



(4.26), by a redefinition of parameters and change of Clifford algebra basis. 



Although, we believe that any general configuration can be generated from 



the triangular solution (4.26) by the change of the Clifford algebra basis or 



unitary rotation, we will present no proof here. 



Space diamond 

If we change the initial setup of the problem, removing the time and promot- 



ing the Hamiltonian (2.1) to the status of Euclidean action, we can repro- 



duce Creutz space diamond fermions [20j from the solution (4.26). Indeed, 
for D' = 4, let us choose the solution with the symmetric "graphene-like" 
location of the Fermi points. 



Kl = k; = -k; = -Kl = 27r/3. 
This choice together with the choice for m/, 

mi = m2 = t, 



(4.30) 



(4.31) 



yields the solution, 

Ti = t. (4.32) 
Using this solution as well as the representation of the D'-dimensional Clifford 



algebra (A. 6) constructed in the appendix, we get a model equivalent to the 



Creutz space diamond proposed in 20 



The resulting Dirac fermion has an internal SU (2) global symmetry due 
to the fact that each Fermi point contributes a full Dirac Fermion rather 
than a Weyl fermion minimally allowed by the Nelsen-Ninomiya theorem 
21 . The pure Dirac fermion can be obtained by an approach similar to 
Kogut-Sussking staggered fermions {22], which, in the case of space dia- 
mond fermions, was successfuUy applied in 23 . The idea is to leave only 



a Weyl fermion at a Fermi point. To do this one leaves only positive chi- 



rality part of the fermionic field in the Hamiltonian (2.2) in odd lattice 
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sites — odd), and only negative chirality part in even lattice sites 

iYli^i ~ even). The chirality in this case should be defined with respect 
to the irreducible Chfford algebra constructed directly from S/. Each Fermi 
point in the low energy hmit contributes a Weyl fermion of different chirahty 
together resulting in a pure Dirac fermion. The same mechanism can be 
apphed also in general context of our model. 



5 Conclusion 

In this work we considered models of fermionic particles on graphs, which 
in the low energy limit flow to Dirac fermion model. A graph satisfying 
this properties is called Dirac lattice. As a necessary condition for this is 
the appearance of linear Fermi points. If the low energy limit admits a 
rotational symmetry, this is also the sufficient condition. In our description 
the graphs are parametrized by the adjacency matrix. Therefore, the Dirac 
lattice condition is expressed as a condition on adjacency matrix. 

We show that the existence of a non degenerate linear Fermi point in D 
dimensional space implies that the energy spectrum of fluctuations around 
it belongs to the space of generators of the D-dimensional Clifford algebra. 

In a special Clifford algebra basis, which we construct, the condition of 
emergence of low energy Dirac fermion is translated into algebraic equations, 
which can be interpreted as polygon closure conditions on the complex plane. 
We provide a general solution to these polygon equations in the case of two 
and three dimensions. The last case is essentially the dimensional reduction 
of the four-dimensional case. We also construct a solution in arbitrary num- 
ber of dimensions, which is superposition of two-dimensional solutions. We 
call this solution holomorphic solution, because it is built using the holomor- 
phic basis in the momentum space. We believe that this solution is also a 
general one. 

It is interesting to point out that the constructions for even and odd 
dimensional spaces are considerably different. This is related to K-theory 



stability of the Fermi points in the two cases 15 . Here we extend the sta 



bility arguments of 18 , allowing us to obtain Dirac particles also in odd 
dimensional spaces, essentially as dimensional reductions of superior Clifford 
algebras. 

One of the objectives of this study was to develop the control of the 
properties of the emergent space-time and simple internal symmetries from 
the fermionic particle on a graph model. Next logical step will be to develop 
the construction of more complicated internal symmetries. In particular, it 
would be interesting to find the criteria for a graph to flow to a physically 
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realistic gauge model like QCD or Standard model. 

As a by-product, the developed teclinique can have some additional ap- 
plications. One important application is lattice simulation. The developed 
approach can be used to construct lattice discretization for fermion contain- 
ing models. Indeed, the constructed graph model is a discretization of Dirac 
fermion and inclusion of gauge or gravity interactions can be reached through 
the deformation of the transition amplitudes. It would be interesting to see, 
however, the computational efficiency of such discretization. 

Another point worth understanding is the relation of our analysis with 
the results of the analysis for crystallographic lattices considered in jl3j . The 
graph described by the adjacency matrix in our case does not necessarily 
satisfy any symmetry and are not restricted to correspond to a physical 
lattice. The physical form of the graph in our case is determined among 
others by the Clifford algebra representation. Finding which representation 
leads to a physical lattice and link them to the crystallographic groups is 
another interesting direction for a future research. 
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A Construction of the Clifford algebra basis 

AU possible Clifford algebra representations are well known [14j. Here we 
show, that for some even dimensional space one can construct a holomorphic 
basis of Clifford algebra generators S/, / = 1, . . . , D' /2 such that all matrices 
of this basis have real entries and satisfy the anti commutation relations 
( pl24l ). 

Since we are paying special attention to the case of two and three spacial 
dimension, let us start with the cases of D' = 2 and, respectively. D' = 4. 
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D' = 2 case 



In the D' = 2 case there is a single nilpotent matrix Si = cr+, where (T+, 
cr^ = O, is the Standard su(2) spin raising operator, 



O 1 
O O 



(A.l) 



It gives rise to the famihar Pauh cr-matrices, which are generators of two- 
dimensional Chfford algebra: 



/3i = -0-1 



(J, 



/^2 = -(72 = i{a^ 



(A.2) 



where (T_ = cr:, 

In the main text, due to Fermi point degeneracy, the Chfford algebra 
generators appear in a reducible representation, 

/3i = /3i ® ^3, /32 = /3i ® I (A.3) 

From these the matrix 7° is introduced according to Equation (3.31), 

7° = P1P2 = ® ^3- (A.4) 

Then, the Dirac gamma matrices are given by muhiphcation of the above 
/5-matrices by 7°, 

7I = -7O/31 = (T2 ® I, 7^ = -7V2 = -fTi ® 0-3. (A.5) 

The general situation here is equivalent to the case of graphene (see an ex- 
cellent review on graphene [8]). The representation (A.5) is reducible. The 



decomposition of matrices (A.5) in terms of a product of an irreducible ma- 
trix and internal symmetry generator can be found in 19 . This situation 
corresponds to minimally doubled Dirac fermion. 



D' = 4 case 



Let us turn to the case of two independent S/ matrices, needed to describe 
the three- and four-dimensional cases. The algebra (3.13) for D' = 4 can be 
realized in terms of two fermionic operators acting on the fermionic oscillator 
basis |si,S2), according to the following rule, 

(A.6a) 
(A.6b) 
(A.6c) 
(A.6d) 



Si 


|0,0) 


= |i,o), 


^2 


|0,0) 


= |o,i), 


si 


|i,o) 


= 0, 


^2 


|i,o) 


= -ll,l) 


si 


|o,i) 




^2 


|o,i) 


= 0, 


si 




= 0, 


^2 




= 0, 
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where Sj = 0, 1 are the occupation numbers. Re-ordering the fermionic 
oscillator basis into the basis e^, i = 1, . . . , 4 as follows, 



we 



e2s2+si+l — S2) , 

represent the matrices S/ in the following form, 



/0 


1 





o\ 























1 


^0 








0/ 



s, = 



/0 





1 


o\ 











-1 














^0 








0/ 



(A.7) 



(A. 



General D' case 

It is not difiicult to generahze the above two constructions to the case of a 
general even dimension D'. In order to do this let us consider the Hilbert 
space of D' /2 fermionic oscillators built from the 'vacuum state' |0), annihi- 
lated by S/, 

|s) = sl^sl^...sljo), S,|0) = 

where the set of indices is ordered as follows, /1 < /2 < ■ ■ 
element for the operator E/ in such a basis is given by, 



(A.9) 
The matrix 



S,|s) 



(A.IO) 



where mg' is the number of filled positions with < /, and s'+/ is the formal 
state which differs from s' by an additional state at position J. If the state 



/ is already filled in s', the right hand side of (A.IO) vanishes automatically. 

By construction, matrices E/ satisfy the algebra (3.13) and have real- 
valued entries. 
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